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Abstract 

Sum rules for the off-shell isovector form factor 7r° — ► 77* are given in terms 
of the pion form factor and the ^yir^ — > n^n experimental data. Similarly, the 
corresponding sum rules for the off-shell isoscalar form factor are given in terms of 
the experimental photon-37r form factor and the 77r — ► 3ir amplitude. For some set 
of parameters, e.g T(p — > ^7) ~ 90KeV, the ir° — > 77 rate and slope parameter of 
the decay tt° — > ^ye + e~ are in good agreement with the experimental data. 
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Dispersion relation is an useful approach to calculate a physical amplitude when non- 
perturbative effects are important. For example, in the presence of a low energy resonance, 
either one has to modify the perturbative results or to use dispersive approach in order 
to satisfy the unitarity relation. The case of the vector pion form form factor, using dis- 
persion relation, was previously discussed M. Sum rules for low energy parameters were 
written down and were found to be in agreement with the experiment data. 

In this article we are interested in calculating the radiative electromagnetic form factor 
of the pseudoscalar meson ir°, for example the form factor of 7r° decaying into 7e + e~ 
or e + e~ — > 7T°7 using dispersion relation. These processes are usually studied by the 
Vector Meson Dominance model @ using the Gell-Mann, Sharp Wagner (GMSW) model 
|| which includes some non perturbative effects. We want to show that the dispersive 
approach enables us to express the calculated results in terms of the measurable quantities 
with a better precision than the GMSW model. Some inacuracies and difficulties of the 
vector meson dominance using the GMSW scheme can thus be avoided ||, [5[ El . 

Our line of approach is similar to the study of the pion form factor The theoret- 

ical calculation was based on the solution of the integral equation of the Muskellishvilli- 
Omnes type || |9|], using the experimental information on the P-wave phase shifts in the 
low energy region up to 1 GeV. This approach yields a r.m.s. radius of the pion too low 
by 10% and the square of the magnitude of the pion form factor at the p resonance (i.e. 
the leptonic width) by 30% |7|. This is also the precision of the VMD model. The 
GMSW model 0, which used the VMD model 0, should therefore be accurate to this 
level. 

On the other hand, if the information on the P-wave phase shifts was combined with 
the experimental magnitude of the time like pion form factor, the r.m.s. of the pion 
could be evaluated by a sum rule which gave a result in agreement with the data to a 
few percents 0. Likewise, if the solution of the Omnes-Muskellishvilli integral equation 
was used with the input of the experimental value of the r.m.s. radius, the accuracy of 
the leptonic width of the p resonance could also be calculated with an accuracy of a few 
percents [|, 0, 0. 

In using dispersion relation for the 3-point function, one has the choice of dispersing or 
taking either one the three particles involved off its mass shell, keeping the other two on 
their mass-shell. For the case of 7r° — > 77* or 7* — > 771"°, one disperses in one of the photon 
mass squared variable. In the following, we denote the dispersed photon as 7* and the 
on shell photon simply as 7. Unlike the GMSW scheme 0, we do not take two photons 
simultaneously off their mass shells, for which dispersion relations are not justified from 
a more fundamental viewpoint. 

Proceeding with the GMSW scheme, one can calculate the ir° — > 77 amplitude which 



was found to agree with the n 77 anomaly |L2[ and also with the experimental data [£L4 



In this calculation, one uses the experimental data on the leptonic widths of p, oj, the 
uj — > 7T7 width and the £77(3) relation relating pir'y vertex to UTC-f. The 737T amplitude 
0, |5], [J in the chiral limit is a factor of 3/2 too large compared with that given by the 
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anomaly JHJ . Although the introduction of a contact term in the 7 — > 37r amplitude could 
solve this problem H, there is a problem with unitarity which will be discussed elsewhere 
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It is therefore desirable to use dispersion relation, in combination with unitarity, to 
solve this problem. 

For the decay process tt° — > 7(^)7* (q), neglecting the electron positron masses, the 
kinematics requires that < q 2 < m 2 . The space-like radiative form factor of ir° can be 
measured by the two photon processes in e + e~ collision, e + e~ — > e + e~ir°, keeping one of 
the exchanged photon almost on its mass shell and the other photon off its mass shell. 
The limit of the space-like momentum transfer is not limited by the tt° mass, but by the 
kinematics of the two photon process. Likewise the time-like radiative form factor of the 
7r°7 can be measured by the process e + e _ — > 7r°7 i.e. the production of 7r° and 7 by 
the single photon exchanged process. There is no kinematical limitation of the time-like 
momentum transfer. We shall study this form factor as a function of the momentum 
transfer squared from —lGeV 2 < q 2 < lGeV 2 . 



1 The Isovector Matrix Element 

The electromagnetic hadronic current is decomposed into two parts, the isovector and 
isoscalar contributions. The lowest number of pion contributions are, respectively, two 
and three pion states. These are the most important states because they are associated 
with, respectively, the p and u vector mesons. Dispersing in the photon mass one must 
take into the contributions of the 2ir and 3tt states with the total angular momentum 
corresponding to the photon or J = states. Because of the G-Parity invariance for 
strong interaction, these two states are orthogonal and we can therefore consider their 
contribution separately. The final result for the form factors 7* —>■ 77r or 7r — > 77* is the 
algebaic sum of their contributions. 

The isovector and isoscalar contributions are denoted, respectively, as F 3 (q 2 ) and 
F°(q 2 ). More explicitly the isovector contribution to the matrix element can be written 
as: 

M 3 (tt( P )) -> 7 (*)7*(g)) = -^ T e,Ke a q T F^{q 2 ) (1) 
and similarly the corresponding amplitude for the isoscalar contribution is denoted by 
F° 7 *(g 2 ). In the following, for convenience, we set q 2 = s. 

At zero momentum transfer s = 0, the off-shell photon becomes the on-shell photon 
and hence the matrix element in Eq. (JJ) is just given by the matrix element of tt° — > 77 
and is given by the chiral anomaly: 

e 2 

where the number of color N c = 3 has been used and f n = 93MeV is the pion decay 
constant and e is the electric charge. It is generally assumed that the matrix element for 
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the physical pseudoscalar 7r° is not very much different from its chiral value, Eq. (|2|). This 
is a plausible assumption, which will be used in this article, because it is well verified for 
7T° — > 27 decay. 

It is clear that Fj? , (s) is an analytic function in the complex s plane. Depending on 
the assumption on its asymptotic behavior and the precision of the calculation, an unsub- 
tracted, once or twice subtracted dispersion relation for FjL* can be written. Similarly its 
first, second ... derivatives at any point except on the positive real axis with s > Am 2 can 
be written in terms of its imaginary part. For simplicity in writing we shall drop, at the 
moment, the superscript and subscript of the funcion and simply write it as F(s). 

Because the form factor F(s) is an analytic function with a cut from Am 2 to 00, the 
n th times subtracted dispersion relation for F(s) reads: 



s n r°° ImF(z)dz 



F(s) = a + a lS + ...a n ^s n - L + — / — , (3) 

7T J Am\ Z n {Z — S — IE) 

where n > and, for our purpose, the series around the origin is considered. Because of 
the real analytic property of F(s), it is real below 4m 2 . By taking the real part of this 
equation, ReF(s) is related to the principal part of the dispersion integral involving the 
ImV(s) apart from the subtraction constants a n . 

The polynomial on the R.H.S. of Eq. (Q) will be referred in the following as the 
subtraction constants and the last term on the R.H.S. as the dispersion integral (DI). The 
evaluation of DI as a funtion of s will be done later. Notice that a n = F n (0)/n\ is the 
coefficient of the Taylor series expansion for F(s), where F n (0) is the nth derivative of 
F(s) evaluated at the origin. 

In theoretical work such as in perturbation theory, the number of subtractions is 
minimal in the sense that it just makes the DI converges. In general, the more subtractions 
we make, the better is the suppression of the high energy contribution which, in general, is 
difficult to calculate. This high energy suppression makes the calculation more reliable, 
but at the same time, makes the calculation less predictive because of the introduction of 
extra parameters. 

ImF(s) can be evaluated using the unitarity relation. Instead of taking a complete 
set of intermediate states in the unitarity relation, we truncate the summation and keep 
only the two pion P state. We shall show below that this approximation is sufficiently 
accurate to analyse the data for s below 1 GeV 2 . 

It is straightforward to show in this case: 

p ( s — 4m 2 W 2 

ImF(s) = — [ - ^_y*( a ) Gl ( a ) (4) 

907T S 1/Z 

where V(s) is the vector pion form factor with the normalisation V(0) = 1 and Gi(s) is 
the P-wave projection of the invariant amplitude for 77r + — > 7r + 7r°: 

M( 1 (k)7T°(q + ) - n + (p + )n°(p )) = ze^q»qlq T _G(s,t,u) (5) 
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where e is the photon polarisation, s = (k + q+) 2 , t = (p+ — q + ) 2 and u = (po — q+) 2 ■ 
Because all particles involved are on shell, one has s + 1 + u = 3m 2 . In the center of mass 
system, in terms of the scattering angle 9, we have: 

t = — - + i(s - ml)^l - 4m 2 /scos0 

St/x^ s ! i 

u = ^ ~(s - m 2 )\Jl - Ami/ s cos 9 (6) 



The partial wave expansion for G(s,t,u) is given as follows [16[ |17|1 : 



G(s,t,u) = J2 G i( s ) P i( cos 0) ( 7 ) 

oddl 



where 9 is the scattering angle and hence 

G 1 (s) = — [ dQ sin 2 9G(s,t,u) 
In terms of the function G(s, t, u) the differential cross section for the process 77r ± 



IS 



da 1 / 2x ( s - 4m 2 ) 3/2 . 2 . , , |2 

(s-m 2 )^ sm 2 g | G(s,t,u) \ 2 (9) 



rfcos^ 1024tt v nJ s 1 / 2 
In Eq. @, the elastic unitarity relation for the vector pion form factor requires that 
the phase of V(s) is the phase of the P-wave 7T7t strong interaction S(s). Similarly, the 
phase of for the P-wave amplitude Gi(s) of the amplitude •yn + — > ty + ti is also S(s). The 
product V(s)*G(s) is therefore real and hence we can rewrite Eq. (f|) as: 

ImF ^ {s) = w, [s ~^ f ' 2 1 v{s) 11 Gl(s) 1 (10) 

where we have restored the subscripts and superscripts. This equation is ambiguous to 
a plus or minus sign. In the following we take the positive sign by convention. Hence 
ImF(s) is only given in terms of the measurable quantities, | V(s) \ and | Gi(s) |. 

There are excellent experimental data on the vector pion form factor but accurate 
experimental data on 77r + — > 7r + 7r° are not available. To proceed, at this stage, one has 
to use some theoretical studies on this process which will be discussed below. 

Besides the 2n, there are higher mass states contributing to the unitarity relation for 
ImF^ t (s). For example, the KK or the An states in the form of the urn or the pirn states 
etc. could become important at an energy above 1 GeV. As long as we are interested 
in the phenomology below this energy region, it may be safe to take only into account of 
the lowest intermediate states which, for the isovector contribution, is the interacting two 
pion state or the p vector meson, and for the isoscalar contribution, the interacting three 
pion states or the uj vector meson. 
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If we want to find out how important is the 2tt contribution in the form of p compared 
with the higher intermediate states, we can write an unsubtracted dispersion relation for 
FjL*(s) and then evaluate this amplitude at s = to compare with the amplitude of 
7T° — ► 27. Needless to say, such an use of an unsubtracted dispersion relation can be 
opened to criticisms because the higher intermediate states are not strongly suppressed; 
we can only hope that, at low energy, their contributions are small compared with the 2n 
and hence can be neglected. This is the spirit of the Vector Meson Dominance model in 
the more precised language of the dispersion theory. 

On the other hand if the e + e~ spectrum in the tt° — > / -/e + e~ process is to be calculated 
the higher intermediate states are then suppressed by a factor s~ 2 (see below) which 
makes the result more reliable. The higher the derivative is at s = 0, the larger is the 
suppression of the high energy contribution. 

We shall parametrize phenomeno logically the pion form factor data below 1.1 GeV by 
the following formula JTT] which agrees well with the experimental data [ IS| , |HS| ]: 



v(] = 1 + 0.16s/ s R 

U 1 - s/s R - ^{(s - Aml)H nn (s) + 2s/3} { ) 

where f n = 0.093GeV, and Sr = 30. 2m\ which gives a good fit to the modulus of the pion 
form factor below 1.1 GeV and the P-wave nir phase shifts in this energy range. H n7T (s) 
is a well-known integral over the phase space factor: 



Is -Ami, V^+Js- Ami s Is - Ami . , 

H^(s) = (2 - L ) + i*y — ^ (12) 

Neither good experimental data on | G(s,t,u) | nor its P-wave amplitude G\(s) are 
available. In fact, the data on the width T(p — > rcy) are in contradiction |20], |22|]. Once 
the experimental data on | Gi(s) \ are available they can directly be used in Eq. (10) to 
calculate the 7r°7 form factor. 

For the time being, it is reasonable to use some theoretical models to parametrize 
Gi(s). This can be done by assuming the elastic unitarity relation and dispersion re- 
lation for this process. A more complicated singular integral equation than that of the 
Muskelishvili-Omnes type is obtained [|lf|. Given the experimental P-wave irn phase 
shifts, its solution, which is related to the T(p — > 717) width, can be obtained by a numer- 
ical method and is a sensitive function of the second derivative of Gi(s) with respect to s 
at s = m 2 . This problem is related to the existence of the contact term in the photon-37r 
coupling 0, ||. Because Gi(s) is an analytic function in the s variable with a unitarity 
right hand cut extending from Am\ to 00, and a left hand cut on the negative s axis. On 
the right cut this function, by unitarity, must have the phase of the strong P-wave nir 
interaction S(s). The left hand cut is due to the t and u channel contributions. From 
general grounds, the P-wave amplitude can be written as a product of the two cuts, the 
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right and left hand cuts. For our purpose, the left hand cut contribution will be, phe- 
nomenologically, represented by the product of a double pole and a term linear in s. We 
have hence: 

q w =v(s)c j^w (l - s/b) (13) 

where c, a, b are real and G® is the chiral anomaly of 7 — > 3n: 

«5 = ^ («) 

Eq. (0), with the real parameters a,b,c, is flexible enough to fit the solution of 
the theoretical calculation. It also gives the correct phase for the P-wave amplitude of 
the process 771- + — > 7r + 7r° which are the P-wave nir phase shifts below 1.1 GeV. This 
parametrization is most accurate in the energy range 0.20GeV 2 < s < 0.8GeV 2 where its 
contribution to the following low energy sum rules is most important. 

Instead of characterizing the results of the theoretical calculation as a function of 
the input second derivative values at s = m 2 , it is more physical to express them as a 
function of the T(p — > tt 1 -/) width which is defined by the value of Eq. (fl3|) compared with 
its corresponding Breit-Wigner formula at s = m 2 . 

In Table 1, we give the values of the parameters as a function of the T(p — > n^) widths. 

Rewriting Eq. ([T0|) in the unit of F 77 as given by Eq. (|[) one has: 

ImF 3 Js) 1 fs-4m 2 W 2 a 2 

; r{ ) = ' 2 ( r I V(s) I 2 c-^(l - s/b) (15) 
F 77 96ttj 2 V s i a + s ) 

Now knowing ImF^, (s) we can write down sum rules and dispersion relation. 

Let us begin with testing the notion of the vector meson dominance, which can be 
translated into the dispersion language, by assuming an unsubtracted dispersion relation 
for F(s). This allows us to find out how much the contribution of ImF^,(s) is in the 
7T° — > 77 amplitude. 

Writing an unsubtracted dispersion relation for F^, (s) and set s = we have: 

c. ( o)=ir^ (16) 

' ' TT J Ami s 

The use of an unsubtracted dispersion relation requires that F^ 7 *(s) — > as s — > 00. It is 
model dependent because the high energy contribution to ImF? n * is only suppressed by 
a factor of s in the sum rule, Eq. (||). 

With this reservation, using ImF^* (s) as given by Eq. (0), with the pion form factor 
I V(s) I given by Eq. (|TTD and with the 'Jtt — > amplitude given in Table 1, the matrix 
element of the isovector tt° — > 77* can be calculated as a function of the T(p — > 7c-f)width. 
The results of our calculation are given on the 5th column of the Table 1. 
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We now want to calculate the energy dependence or the pair e + e~ in 7r° — ► / -/e + e~ at 
low energy. This can sufficiently be done by evaluating the first derivative of FjL»(s) at 
the origin: 

1 ' 7T J 4ml S 

and if more precision is needed one could include higher derivatives. For example the 
second derivative is given by: 



7T J 4ml S 3 



The sum rule for the first derivative at s = of FA,* is quite reliable because the unknown 
high energy contribution is suppressed by the factor 1/s 2 in Eq. (|T7j). The slope is usually 
expressed in terms of of the slope parameter j3 l defined as /3*s/m° 2 with i = 0, 3 for, 
respectively, the isoscalar and isovector contributions. 

In table 1, the slope parameter for the isovector contribution is also given in the 5th 
column. The precision of the experimental data do not allow us to compare the second 
derivative of *. We shall compare the values of the matrix elements and slopes after 
our calculation of the isoscalar contribution. 

In evaluating the RHS of Eq. flU) or Eq flT7|), the P-wave phase space factor in Eq. 
(|I0|) suppresses strongly the low energy contribution of the modulus of the amplitude 
77r — ► 7r + 7r~. They are therefore very insensitive to the low energy behavior of this 
amplitude. The well-known problem of reconciling Vector Meson Dominance and the 737r 
chiral anomaly is not a problem here. It will be discussed in a separate study 1T5 . 



2 The Isoscalar Matrix Element 

Dispersing in the isoscalar photon mass, we have to take into account of the lowest inter- 
mediate state in the unitarity relation, the 3tt states. The total isospin of the 3 pions is 
1 = and its total angular momentum is 1~ hence they are spacely completely antisym- 
metric. The elastic unitarity contribution to the isoscalar amplitude is: F^*(s). 

ImF° f 0) =< 7 | 3tt >< 3vr I 7 7T > (19) 

Again, after integrating over the internal variables for the 3 pion system, we can establish 
the phase theorem for < 7 | 3tt > as in the case of the 2 pion contribution to the isovector 
form factor. The RHS of Eq. (|19D is therfore real and we can rewrite it as: 

ImF° r (s) =|< 7 I 3vr >||< 3tt | 7 7r >| (20) 

Hence ImF® * (s), similar to the isovector case, can be directly obtained from experimental 
data. There are good data for the process 7* — > 3n which shows a complete dominance 
of the lj contribution. 
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The process 77r — > 3n has not been measured. However, the 3 pion resonance is such 
a narrow one that it is possible to approximate it accurately by the u meson. This is also 
true for the photon 3 pion vertex. 

Defining the photon u coupling as em^/g^ @ with the photon invariant mass squared 
equal to the to mass the g u can be determined from its leptonic width: Y{u — > e + e~) = 
(Ana 2 /3)(m UJ /gl). The u -> tt°7 width is given by: T(lu -> tt°7) = g^E? ( /(l2ir). Using 
the experimental values, m w = 782MeV, T w = 8.41 ± .09MeV and the branching ratio 
r^Try/r^ = 8.5 ± .5%, one can easily verify: 

F° 7 = 0.50 ± 0.02F 77 (1 + -^-) (21) 

hence the a; contribution is just half of the 7r°77 anomaly, Eq. (0). The imaginary part 
of the isoscalar form factor can be straightforwardly deduced from this result: 

where T(s) is strongly dependent on s due to the three pion phase space. We shall neglect 
altogether the real part of the u self energy correction. 



3 Sum of the isovector and isoscalar matrix elements 

Knowing the imaginary part of the form factor, Eqs. (KB, ^) and using dispersion relation, 
there is no difficulty in calculating the form factor in the interval —lGeV 2 < s < lGeV 2 . 
The energy dependence of the n° — > 77*, as will be shown below, requires us to take the 
same sign for the imaginary part of the isoscalar and isovector contributions. 



3.1 7T° — > 77 and 7r c 



77 



Let us first examine the 7r° — > 77 and 7r° — > 77* amplitudes. The observable matrix 
element _F 77 *(s) is the sum of the isoscalar and the isovector matrix elements. They 
cannot experimentally be separated from each other. We want to find out how much 
the two and three pion intermediate states (in the form of u and p) contributing to the 



tt" — > 77 amplitude. Using Eq. (|2l ) and the results given on the fifth column of the 
Table 1, the total matrix elements are given in Table 2. They should be compared with 



experimental value [14, 21 



F 77 ,(s) = 1 + (0.032 ± .004)s/m^o 



(23) 



Comparing this result with those given by Table 2 it is clear that the matrix element 
corresponding to T(p —>■ ir / -f) = 92KeV is favored. The experimental errors are, however, 
large. 
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Because we have more faith in the slope sum rule than the unsubtracted dispersion 
relation due to the suppression of the high energy contribution the prediction on the 
r(p — > TT'-f) width can reliably be made once the slope is better measured. 



3.2 Time-Like and Space-Like 7r°7 Form Factor 

Using the imaginary part of the isovector and isoscalar contributions to the form 
factor F 77 *(s) as given by Eqs. (]T5|, we can calculate the s-dependence of the form 



factor from —lGeV 2 < s < lGeV 2 with reliability by the following subtracted dispersion 
relation: 

iV(s) _ 1 + s r°° ImF ir {z)dz 



F 77 7i J Ami z(z — s — ze) 



where JmF Tr «(s) is given by the sum of the RHS of Eqs. flloT , y^A)- 

Because of the overwhelmed dominance of the uj contribution to the form factor F 77 (s) 
around the uj — p region, there is little hope to get physics out of the p and uj region 
unless the accuracy of the cross section measurement can be improved to a few percents. 
Likewise, an improvement in the determination of the space like behavior of the form 
factor will enable us also to determine accurately the pTcy width. 



4 Conclusion and Acknowledgements 

The main point of this paper is to show the usefulness of the dispersion relation approach 
which has not effectively been used in the literature for this problem. It allows us to 
relate various experimental measurements without making strong assumptions. The non 
perturbative effects due to strong interaction are straightforwardly handled by dispersion 
relation and unitarity. 

Part of this is done at the Institute for Nuclear and Particle Studies at the KEK High 
Energy Accelerator Research Organization. The author would like to thank Professors 
M. Kobayashi and H. Sugawara for hospitality. 
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r(p -> 7T7) 


c 


a 


b 


Isovector Matrix Element (7r — > 77*) 


60 


1.01 


80 


80 


0.40 + 0.0145s/m2 o 


77 


0.97 


80 


100 


0.44 + 0.0154s/m2 


92 


1.06 


50 


oo 


0.49 + 0.0171s/m2 


110 


1.04 


60 


oo 


0.55 + 0.0193s/m2 () 



Table 1: Parameters describing the P-wave amplitude for 771 — > txtx and the isovector 
matrix element in the unit of the anomaly _F 77 Eq. (0); a and b are in the units of 

= 1. r(p — ► 7T7) are in the unit of KeV. 



r(p 7T7) 


Total Matrix Element (ir — > 77*) 


60 


(0.90 ± .036) + (0.0292 ± .0012)s/m2 


77 


(0.94 ± .037) + (0.030 ± .0012)s/m2 


92 


(0.99 ± .04) + (0.032 ± .0013)s/7t4 


110 


(1.05 ± .04) + (0.034 ± .0014)s/m2 



Table 2: Total matrix element is the sum of isovector and isoscalar matrix elements in 
the unit of the anomaly F 77 Eq. (fj). T(p — > 717) are in the unit of KeV. The errors come 

from those of the isoscalar contribution. 



10 



References 

[1] T. N. Truong, preprint (|hep-ph/000127l| ). Talk given at the 5th workshop in QCD, 
Villefranche-sur-Mer January 2-7, 2000; Proceedings (to be published). 

J. J. Sakurai, Ann. Phys. (N.Y) 11, 1 (1960); M. Gell-Mann and F. Zachariasen, 
Phys. Rev. 124, 953 (1961); Y. Nambu and J. J. Sakurai, Phys. Rev. Lett. 8, 79 
(1962); 8, 191(E) (1962). 

M. Gell-Mann, D. Sharp and W. G. Wagner, Phys. Rev. Lett. 8 261 (1962). 
S. Rudaz, Phys. Rev. D 10, 3857 (1974); Phys. Lett. 145B, 281 (1984). 
O. Kaymakcalan, S. Rajeev and J. Schecter, Phys. Rev. D 30, 594 (1984). 
M. Bando, T. Kugo and K. Yamawaki, Phys. Rep.164, 217 (1988). 
T. N. Truong, Phys. Rev. Lett. 61, 2526 (1988). 

N. I. Muskhelishvili, Singular Integral Equations (Noordhoof, Groningen, 1953). 
R. Omnes, Nuovo Cimento 8 316 (1958). 
T. N. Truong, preprint ( hep-ph/0101345 ). 
Le viet Dung and T. N. Truong, Report No [hep-ph/9607378. 



J. S. Bell and R. Jackiw, Nuovo Cimento 60, 47 (1969); S. L. Adler, Phys. Rev. 177, 
246 (1969). 

S. L. Adler, B. W. Lee, S. B. Treiman and A. Zee, Phys. Rev. D 4, 3497 (1971); M. 
V. Terent'ev, Phys. Lett. 38B, 419 (1972); J. Wess and B. Zumino, Phys.Lett. 37B, 
95 (1971); R. Aviv and A. Zee, Phys. Rev. D 5, 2372 (1972). 

Particle Data Group, Eur. Phys. J. C 3, 1 (1998). 

T. N. Truong (under preparation). 

M. Gourdin and A. Martin, Nuovo Cimento 6, 1102 (1960); H. S. Wong, Phys. Rev. 
Lett. 5, 70 (1960). 

T. Hannah (to be published). 

L. M. Barkov et al. Nucl. Phys. B256, 365 (1985). 

ALEPH Collaboration, R. Barate et al, Z. Phys. C 76, 15 (1997). 



11 



[20] L. Capraro et. al, Nucl. Phys. B288, 659 (1987); J. Huston et. al. Phys. Rev. D 33, 
3199 (1986). 

[21] H. J. Behrend et. al, Z. Phys. C 49 401 (1991). 
[22] S. I. Dolinsky, et. al. Z. Phys. C 42 511 (1989). 



12 



